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We developed analytical technique for extracting the phase, visibility and am- 
plitude information as needed for interferometric astrometry with the Space 
Interferometry Mission (SIM). Our model accounts for a number of physical 
and instrumental effects, and is valid for a general case of bandpass filter. We 
were able to obtain general solution for polychromatic phasors and address 
properties of unbiased fringe estimators in the presence of noise. For demon- 
stration purposes we studied the case of rectangular bandpass filter with two 
different methods of optical path difference (OPD) modulation - stepping and 
ramping OPD modulations. A number of areas of further studies relevant to 
instrument design and simulations are outlined and discussed. 
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Introduction 

SIM is designed as a space-based 10-m baseline Michelson optical interferometer oper- 
ating in the visible waveband (see Ref. 1 for more details). This mission will open up 
many areas of astrophysics, via astrometry with unprecedented accuracy. Thus, over 
a narrow field of view SIM is expected to achieve mission accuracy of 1 //as. In this 



mode SIM will search for planetary companions to nearby stars by detecting the as- 
trometric "wobble" relative to a nearby (< 1°) reference star. In its wide-angle mode, 
SIM will be capable to provide a 4 fias precision absolute position measurements of 
stars, with parallaxes to comparable accuracy, at the end of a 5-year mission. The 
expected proper motion accuracy is around 3 /xas/yr, corresponding to a transverse 
velocity of 10 m/s at a distance of 1 kpc. 

The SIM instrument does not directly measure the angular separation between 
stars, but the projection of each star direction vector onto the interferometer baseline 
by measuring the pathlength delay of starlight as it passes through the two arms 
of the interferometer. The delay measurement is made by a combination of internal 
metrology measurements to determine the distance the starlight travels through each 
arm, external metrology measurements that determine the length and local orienta- 
tion of the baseline, and a measurement of the central white light fringe to determine 
the point of equal optical pathlength (see Figure QJ). 

The primary motivation for the work was the idea to use the averaged and bias- 
corrected complex phasors to estimate the external optical pathlength difference for 
the incoming polychromatic light. This approach is applicable to the science interfer- 
ometer phase measurement, for which it is not necessary to remove the bias from each 
estimate, just as long as it is removed from the averaged estimate. This method was 
applied in Ref. 121 to obtain unbiased fringe visibility estimates. Under certain condi- 
tions (dc intensity values and visibilities approximately constant during integration 
time) this method allows to use the average values of the phasors for estimating the 
average phase (see discussion in Refs. EHU). A definition for the complex visibility 
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Fig. 1. Basic geometry of light propagation in the stellar interferometry. 

phasors is stemming from the form of a complex visibility function, V = ve^, where 
v is the visibility and is its phase. Decomposing this expression onto real and imag- 
inary parts as V = X + jY, one obtains the complex visibility phasors, X = v cos <fi 
and Y = t>sin0. SIM will be able to effectively determine both visibility and phase 
of the fringe, but for the astrometric purposes the phase must be determined to a 
much higher accuracy. (The most stringent SIM requirement in this regard is the 
average phase error over 30 s integration time corresponds to a path-length error of 
approximately 30 pm.) Phase determination in the presence of noise is a non-linear 
estimation process. Even in the monochromatic case it requires careful approach to 
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averaging and correcting for biases in the data. Most of the algorithms for phase- 
shifting interferometry are designed for monochromatic light, and there is typically 
a match between the stroke of the modulating element and the wavelength of the 
light (see, for example, Refs. El for a discussion of various algorithms). However, be- 
cause SIM uses a dispersed fringe technique, this match cannot be maintained over 
multiple channels. A number of modifications to the existing four-bucket algorithms 
were designed to address specific problems relating to Palomar Testbed Interferome- 
ter (PTI)P and the Keck Interferometer. Also, there is nothing inherent that requires 
the match between wavelength and stroke, nevertheless most analysis has been based 
on this assumption. It was shown in Ref. HJthat at low light levels, most of these algo- 
rithms become biased when there is a mismatch. The magnitude of this bias is quite 
significant with respect to the instrument requirements, thus making these algorithms 
unacceptable for SIM. 

This paper discusses analytic model developed for the white light fringe data 
extraction. Our goal here is to establish functional dependency of the white light 
fringe parameters on the properties of incoming light as well as the instrumental 
input parameters. We will show that this approach is applicable for establishing the 
unbiased estimators for the case with low light levels. This method enables on to 
analytically analyze the noise propagation properties. The importance of this feature 
comes from the fact that integration time on stellar targets accounts for a significant 
portion of the mission time and can be reduced by the use of processing techniques 
with reduced error variance (see Refs ITHTUI for more details). The analytical from of the 
fringe parameters may be helpful in studying different properties of the instruments, 
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especially their contribution to the accuracy of astrometric delay measured by SIM 
science interferometer.^^ 

The problem of interference of electromagnetic radiation is well studied and ex- 
tensive number of publications on this subject (specifically related to stellar interfer- 
ometry) are available (see Refs. llOp^HTvll and references therein). However, because of 
complexity of this problem in a general case of polychromatic light, most of the current 
research is done numerically. While numerical studies have proven to be extremely 
valuable in analyzing the interference patterns and are very useful in addressing var- 
ious instrumental effects, the analytical methods provide the much needed critical 
understanding of the white light interference phenomena. It will be demonstrated be- 
low that analytic solution may be used as a tool to study the complex interferometric 
phenomena on a principally different qualitative levelP 

In this paper we derive analytic model that may be used to describe photo- 
electron detection process. We analytically describe the physical and instrumental 
processes that are important in estimating the fringe parameters (i.e. intensity of 
incoming radiation, its visibility and the phase of the fringe). Effects that are not 
included in the model are due to polarization of both incoming light and the instru- 
mental throughput, effect of the wavefront-tilt, low frequency vibrations, drifts, jitter, 
etc. Consideration the size of this paper, we plan to address these issues elsewhere. 

The paper is organized as follows: In Section^ we review the description of inter- 
ferometric pattern in the case of monochromatic radiation. In Section |2] we develop a 
model for the interferometric pattern registered by a CCD detector in the polychro- 
matic case. Our model accounts for the effects of the instrumental throughput, beam 
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splitter and quantum efficiency of the CCD. We also discuss the spectral channels 
with narrow bands designed to filter the polychromatic light. In Section El we intro- 
duce parameterization for the polychromatic fringe pattern and define the quantities 
that are forming the astrometric signal on the CCD. Specifically, we derive solution 
for the white light fringe equation in the general case. In Section H] we present general 
analytic solution for complex visibility phasors and will discuss a noise suppression 
approach. In Section |3] we develop technique for studying the case of a rectangu- 
lar bandpass filter. We also obtain functional dependency of our solution in the two 
cases of OPD modulation, namely the stepping and ramping modulations. In Section 
El we present conclusions and recommendations for future studies of accurate fringe 
reconstruction. In order to make access of the basic results of this paper easier, we 
will present some important calculations in the Appendices. Thus, in Appendix 1X1 we 
discuss two possible definitions for the fringe phase and justify the choice we made in 
the paper. In Appendix El we develop approximation for the complex fringe envelope 
function. In Appendix we present a general solution for the instrumental contribu- 
tion affecting the fringe parameters in the case of the rectangular bandpass filter and 
stepping OPD modulation. 

1. Monochromatic Fringe Pattern 

The problem of interference of a monochromatic radiation presently is well studied 
(see ReffT5l and references therein). Here we would like to review information that 
will be necessary for discussion of the fringe phase extraction process. 

We assume that a plane electromagnetic wave, E, which is coming from infinity 
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simultaneously on the two arms of an interferometer, has the following form: 

E = E e j (*"-*■*) , (1) 

with Eq is a constant vector. We will be using a nomenclature where a wavenumber k 
relates to the wavelength as follows k = After passing through the interferometer, 
the two beams E\ and E 2 , are combined at the detector to produce interferomet- 
ric pattern (see Figure . The corresponding fringe pattern is due to the coherent 
addition of the two light-beams, E — E\ + E 2 , and it may be expressed as follows: 

1 = {E ■ E*) t]me = E 2 m + E 2 02 + 2(E 1 E 2 ), (2) 

where * denotes a complex conjugate quantity and (..) t ime denotes a time-averaged 
quantity. The first two terms on the right-hand side of this equation are constant 
intensities of light in the two beams, X\ = Eq X and T 2 = Eq 2 . The third one is 
the interferometric pattern which, in ideal situation, depends only on the optical 
path difference between the two beams and the wavelength of the radiation, namely 



— 02 = -r^o = kdo. The resulting intensity of radiation on is given 



l(k) = 2b (1 + Vcoskd ), (3) 

where we denote the constant part in the intensity pattern as X = I\ + X 2 . Also, 
the quantity V is the visibility of the incoming light, which, in the case of the beams 
intensities mismatch X\ ^ X 2 , is given b} * 13 * 16 * 



= 2y/ZA v = 2\E m \\E Q2 \ v 
Xi+X 2 El 1 + El 2 



The factor V (0 < V < 1) in this equation is the true source visibility (or the 
fringe contrast). In the case of monochromatic radiation from an un-resolved source, 



when the intensities in the two arms of the interferometer are of equal amplitudes, 
or \Eoi\ = l-Ebil = Eo, the constant intensity becomes X\ = I2 = \T$. Therefore, 
visibility reduces simply to one, V = 1. 

A. Fringe Modulation 

As seen in Figure El monochromatic pattern is a simple one. In practice, the studied 
source may be interfere metrically resolved, thus leading to a reduced visibility V ^ 1. 
Moreover, the visibility V, the constant phase difference <po and the constant intensity 
Iq are generally not known. To find them one modulates the phase difference with 
some known function, say x(t), as below: 



where effect of the beam splitter (which is true for a Michelson stellar interferometer) 



The primary goal here is to determine the phase 0o by modulating the internal 
optical path difference x(t). By doing so, one finds the exact value of internal delay 
^int = x(t') that would exactly compensate the initial offset or k(d — x- mt ) = 0. Then, 
having determined the external delay do = x- mt one can determine the source position 
from equation do = (b ■ s), where b is the interferometer baseline vector and s is the 
source position on the sky. In practice, this is done in a global astrometric solutions 
discussed in details in Refs lTTI IT^l and not addressed here. 

As shown in the Figures El and 0] there may be different ways of modulating 
the internal OPD in an interferometer. In particular, Figure El demonstrates the case 
when the OPD is modulated linearly. Thus, the upper plot shows a typical phase 




(5) 



brings additional | phase shift (discussed in Section I2b|) ISMSiai 



8 



Optical pathlength difference, x 

Fig. 2. Monochromatic fringe. Note the visibility is set V = 1 and the initial 
phase offset chosen to be 0o = f . 

change with ramping over 8 equal temporal bins with duration of 10 ms each. The 
lower plot shows the same monochromatic fringe shown in Figure El as a function of 
OPD that is modulated by ramping over a wavelength. Note the difference between 
this case and the case shown in Figure El where the upper plot shows a typical phase 
modulation provided by the stepping OPD modulation with the stroke stepping over 
8 equal temporal bins with duration of 10 ms each. The lower plot in this Figure 
demonstrates behavior of monochromatic fringe as a function of OPD with modulated 
stroke stepping in equal steps over a wavelength as shown in the upper plot. 

The photon count on the detector, N, is proportional to the intensity, N oc X. 
Therefore, by collecting photons Ni(k) ~ T(k,t{), coming at the detector at a certain 
time intervals £j G [t~,tf]), i G [1, .., N] (or by integrating Eq. (JHJ) over dt from £j_x 
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Fig. 3. Upper plot shows a typical OPD modulation stroke with ramping over 
8 equal temporal bins with duration of 10 ms each. Lower plot is a monochro- 
matic fringe as a function of OPD that is modulated by ramping over a wave- 
length (initial phase offset is 0o = f - same as in Figure |2J). 
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Optical pathlength difference, [bins] 



Fig. 4. Upper plot shows typical OPD modulation stroke with stepping over 
8 equal temporal bins with duration of 10 ms. Lower plot is a monochromatic 
fringe as a function of OPD modulated in 8 equal steps over a wavelength, as 
shown in the upper plot. (The fringe parameters the same as in Figure El) 
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to tj), one forms the system of equations to determine the unknown quantities Iq, V 
and 0o- These time intervals correspond to different values of OPD, x(t) E jxj, i E 
[1, .., N]\, therefore, observational equation in the case of monochromatic light and 
ramping OPD modulation takes the following form: 



where sinc(-)-function is given as usual sine z = sin z/z; v is the constant velocity of 
OPD modulation stroke, x(t) = vt and At; = t\ — is the integration time for the 
z-th temporal bin. 1 

Eq. (JUJ) is the most studied equation when estimating the fringe parameters in 
the monochromatic light approximation. Solution to this equation is quite straight- 
forward and it was extensively discussed in literature (see, for example, Ref|2J[T5J). 
Usually, Eq. © is represented in a matrix form as Zj = A ia X a , where indexes i and 
a running as i 6 {1,...,N} and a E {1,2,3}. Vector X a is the to-be-determined 
phasors vector given as X a = (Z ; IqV cos^ Z o ^ sm 0o) ■ Matrix A = A ia = 
(l; sine [~kv At;] sin kx^ sine [\kv Arj] cos kxij is the 3 x ]\f matrix of 3D rotation 
in the phase space. A solution to this equation is given by X a = A^ ta Ii, where 
At = (A T A) _1 A T , with A^ being the pseudo-inverse of A. This set of equations may 
be solved uniquely only in the case when N = 3. For all other cases, when N > 3, 
the obtained system of equations is over-determined, and one obtains a least-squares 

solution by constructing pseudo- inverse matrix A^PEEHWe will discuss an optimally- 

1 Note, by taking the limit v — > in Eq. © (i.e. sivc[^kv Arj] — > 1) one recovers the case of 



stepping OPD modulation with a familiar simple form of observational equation: Ii{k) = + 
Vsin (0o + kxi) J . 
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weighted, noise optimized solution to this set of equations in more details in Section 
while dealing with a more complicated case of the polychromatic light. 
When dealing with a polychromatic light content in the wide bandwidth, one 
either i) forms a light with a narrow spectral width, such that effects of polychromacity 
(discussed further) are negligible, or ii) disperses the light beam on a large number 
of spectral channels, such that monochromatic approximation is valid within each 
channel.^ As a result, most of the current algorithms and simulations for stellar 
optical interferometry are based on the properties of the monochromatic light. This 
is a good approximation for some of existing testbed configurations that use as many 
as 80 spectral channels for dispersed light. Nominally the SIM flight system will use 
four to eight channels for guide interferometers. Because of the large bandwidth of 
each channel (87.5 nm), the quasi- monochromatic assumptions are not valid, and 
modifications to the algorithms are necessaryP^ In the following Section we will 
introduce a method designed to address this issue. 

2. Modeling Observables for a Polychromatic Fringe 

The observational conditions in the case of polychromatic light are significantly al- 
tered compare to the simplicity of the monochromatic situation discussed in Section 
^ Thus, Figure 03 shows a general behavior of harmonic signals with a different fre- 
quencies. The left plot in the figure shows three sinusoidal monochromatic waves with 
the same initial phase offset 0o — tt/5. Depending on the wavenumber, each wave, 
as a function of OPD, produces different fringe pattern. Note that a combination of 
these three waves (a simple model of a polychromatic fringe) also produces a har- 
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monic signal, but its shape is drastically different from the initial one (shown on the 
right plot). For the expected SIM finite bandwidth the wavenumbers of the interfering 
light may be different as much as twice from each other (i.e. the SIM wavenumber 
bandwidth is k G [450,950] nm). Thus, in general, only at zero OPD (i.e. for the 
white light fringe) the interferometric pattern would have maximum intensity. 

In this Section we will derive an equation to describe the interference pattern of 
a polychromatic fringe as needed for the SIM CCD detector read-out. 

A. Modeling Observables in the Wide Bandwidth Case 

As was mentioned earlier, the case of describing interferometric pattern that involves 
a finite bandwidth - is a more complicated one. In general all the quantities involved 
are complicated functions of a wavenumber. Note that, contrary to the Eq. (0), not 
only the phase of the electromagnetic wave depend on time and wavelength, but the 
same true on it's amplitude 

E(k,t) = E Q (k,t)e j{ut - kx \ (7) 

with E being the wave's amplitude vector that depends both on the wavelength and 
time. This expression shows a wave-packet that describes individual contribution to 
the polychromatic light that is due to a particular monochromatic constituent of it. 

A way to describe this process is to collect all constituents of polychromatic light 
at different wavelengths that are present in the incoming electromagnetic wave. It 
is convenient to express the intensity of the polychromatic radiation as a coherent 
addition of the individual wave-packets ([TJ with a different frequencies.^ Denoting 
Eq\ and Eq2, to represent the light coming onto a detector from the two arms of 
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Optical pathlength difference, x Optical pathlength difference, x 

Fig. 5. Left plot shows three independent monochromatic waves presented as 
a functions of the OPD. The phase of the fringes is changing as = 0o + kiX, 
with external delay corresponding to a phase difference of 0o — f an d three 
wavenumbers chosen ki = ko, k<i = 1.5ko, k% = 2ko (thus the bandwidth is 
Ak = ko), where k$ is a reference wavenumber. The right plot represents the 
interference pattern of polychromatic light composed from the same three si- 
nusoidal waves shown on the left plot. Note the drastic change in the character 
of the interferometric pattern. 
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interferometer, this procedure may functionally be expressed in the following form: 
J(M) ~ ((E 01 {k,t)+E 02 {k,t)) ■ (E* 01 (k,t) + E* 2 (k,t))) = 

\ I time 

= 4\(A;) + E 2 2 (A;) + 2|E 01 (A;)||E 02 (A;)|7(A;)cos(<l>(A;)+A^ 12 (fc,t)), (8) 

where A0 12 (fc, t) is the external phase difference between the two arms of interferom- 
eter and ^f(k) = V(k)T(k) = 7(A;)e 3 '*'*) is the complex coherency factor. 

This complex coherency factor includes both - the true complex visibility of the 
source, V(k), and the complex instrumental transfer function denoted as T(k). The 
complex source visibility function is given in its usual form 

V(jfc) = V(fc)e^ v(fc) , (9) 

with V(k) being the amplitude and <pv(k) the phase of the true source visibility. 
According to Van-Cittert-Zernike theorem, namely this quantity is connected to the 
true radiation emitted by a star. Specifically, this theorem allows to calculate mutual 
intensity observed on a surface some distance from the source.^ The instrumental 
complex transfer function, T(k), may also be presented in a similar manner 

f(k)=T(k)e j,t,T{k \ (10) 

with T(k) being the amplitude of this function and </>r(&) its phase. Therefore, the 
complex coherency factor j(k), that is present in the equation (JHJ), has the following 
functional dependency: 

j(k) = V(k)f(k) = V(k)T(k)e j[Mk)+Mk)] = j(k)e jm , (11) 

where combined apparent visibility and phase are given by the expressions j(k) = 
V(k)T(k) and $(&) = <f>v(k) + (pr(k) correspondingly. 

16 



The modulated phase difference Aip^k, t) in the Eq.(jHJ) may be expressed in 
terms of the internal pathlength difference: 

A<p 12 (k,t) = kx(t), (12) 

with x(t) = £i{t) — 0,2{t) being the internal pathlength difference between the two 
arms of interferometer (see Figure Q). 

Let us define individual spectral densities of photon flux in each arm of interfer- 
ometer as X\ ~ E^ik) and X 2 ~ EQ 2 (k). This allows us to present the intensity of 
incoming radiation at a particular wavenumber k as 

I(k,t) ~ Ji(fc) + J 2 (fc) + 2 y /l 1 (k)l 2 (k)i(k) cos ($(*) + kx(t)). (13) 

It is naturally to define the total spectral density of photon flux of light approaching 
the detector as X (k) = Xi(k)+X 2 (k). Note that individual spectral densities of photon 
flux in the two arms of interferometer, X\ and Z2, may be different. To account for 
such a mismatch we introduce apparent visibility V(k) 

*M = t^!t^ w = nw™. (14) 

Xiik) +X 2 (fc) 

The resulted expression for the spectral density of photon flux of the incoming radi- 
ation at a particular wavenumber k takes the following form: 

X(k,t) =X Q (k)(l + V(k)cos($(k) + kx(t))\ (15) 

Furthermore, the energy density of incoming polychromatic radiation per a unit 
wavenumber may be obtained by integrating the photon flux at the detector Eq. (|15j) 
over the wavenumber space dk: 

dX(k,t) = X (k)(l + V(k)cos(®(k) + kx(t))^jdk. (16) 
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Expression (|16[) defines the energy density of the interfering light per a unit wavenum- 
ber. The quantity that is proportional to it, namely the photo-electron count, will be 
registered by a CCD detector via a photo-electron emission process and will be dis- 
cussed in Section I2CI 

B. Effect of a Beam Splitter 

To complete the formulation of our model we need to account for yet one more 
element that is of crucial importance for a Michelson stellar interferometer - the 
beam splitter!^ It is well-known that one of the figures of merit when considering 
various beam splitter designs for the astrometric beam combiner, is the variation of 
phase as a function of wavelength between the interfering beams at zero OPD. Thus, 
ideally, one would like zero phase difference for all wavelengths at zero OPD since 
this would produce a fringe maximum simultaneously at all wavelengths (or fringe 
minimum for a Michelson stellar interferometer). Another desirable case is one in 
which the phase difference varies linearly with frequency. In this case one can still 
achieve a simultaneous fringe maximum at all wavelengths; however, for this case 
there is a constant correction that must be applied to the calculated external delayP 
This effect is important because it may produce a significant contribution to the 
OPD. We may include this effect in a most general way — the wavenumber dependent 
function f(k) added to the desirable 7r/2 effect of an ideal beam splitter. Thus, the 
total effect of the beam splitter may be modeled as 

= £ + /(*), (17) 
where f(k) is a slow varying function of a wavenumber. This function may be ap- 
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proximated up to the second order around some central frequency of bandpass (the 
whole 80 channels pass), k c , as 

f(k) = f(k c ) + f'(k c )(k - k c ) + \f"{k c ){k - k c f + 0(Ak 3 c ). (18) 

It may be showrPthat function f(k c ) will have a direct impact on the phase accuracy 
estimation by shifting the phase by a constant value. In addition, the function f'(k c ) 
will change the envelope function and, in general, will produce a non-linear contribu- 
tion to the phase. f"(k c ) will have impact on both - linear contribution to the phase 
and the non-linear one; first it comes as a correction to the envelope function and 
then to the phase. Even though, the magnitude of the effects of /" is smallest among 
all, the parameters /, /', /" are all of importance. The corresponding effects may be 
well modeled, depending on the optical properties of the beam splitter, but the final 
answer would come probably form calibration. 

Finally, without loosing generality, the effect of beam splitter may be accounted 
for as an additional phase shift to the argument besides the usual 5§b s (k) = ~ and, 
thus leading to a new definition for the fringe phase <p{k) in the form 

_ ^ k ) + s^ bs (k) = ^k) + ^ + f(k) = ^ + 4>(k)- (19) 
Hence, the light intensity Eq. (JTfij) may be given as follows: 

dX{k,t) =l (k)(l + V(k) sin ((j)(k) + kx(t))^jdk, (20) 
where we accounted for the nominal | phase shift due to the beam splitter. 
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C. Detector's Photo-electron Counts 

A CCD detector is responding a quantity that is closely related to the intensity of 
radiation, namely the incoming energy which is given as d£(k,t) = dX(k,t)dt, or 

d£(k, t) = F Q (k)Ta(k) (l + V(k) sin (<f>(k) + kx{t)))dkdt, (21) 

where Fo{k) is a dimensionless factor representing the total instrumental through- 
put Usually a detector is optimized to work at a certain wavelengths better then 
at the others. This property may be qualitatively described by, so called, the quan- 
tum efficiency of a detector, a(k). Thus, the quantum efficiency of a CCD detector is 
conventionally defined as J\f(k,t) = a(k)£(k,t), therefore the density of the emitted 
photo-electrons per a unit area may be presented by the following expression: 

dAf(k,t) = a(k)F (k)l (k) (l +V(k) sin (<f)(k) + kx{t))^jdkdt. (22) 

Note that quantum efficiency of the detector is a function of a wavelength a(k). This 
dependency will not addressed in the present study, but will be explored elsewhere. 

The total number of photo-electron counts, N, depends on the collective area of 
the detector in accord dN(k,t) = Af(k,t)dA. Note that the total photon count is a 
function of power of radiation approaching the detector which is given as dV(k,t) = 
S(k,t)dA. Therefore, one obtains 

dN{k,t) = a(k)J r (k)l (k) (l + V(k) sin (<f>(k) + kx(t))^j(n ■ r A )dAdkdt, (23) 

where n is the direction of the wave falling on the detector and dA = dAfX is the vec- 
tor of the collective area of the detector. Integration of this equation over the collective 
area depends on the properties of the experimental setup, notably on the orientation 
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of the collective area with respect to incoming light. Thus, a small geometric mis- 
alignment in the system may be responsible for the uneven illumination of different 
pixels on the detector. Additionally, the wavefront tilt may produce a measurable 
contribution to the fringe parameters estimated with the help of Eq. ([23)1 . 

In SectionElwe will introduce a concept of a filtered polychromatic light for which 
notation dA# will be designated for the collective area of the £-th spectral channel 
(or area of illumination on the detector designated for a particular spectral channel). 
Assuming that this area is small, and neglecting divergence of the radiation in the 
instrument, we may integrate Eq. ()23|) over dAg. Note, that this operation, together 
with the quantum efficiency of the detector, mathematically may be "folded" into 
the definition for the bandpass filter (also discussed Section EJ). Therefore, we can 
designate to the filter not only the function of control over the allowed band-pass 
of incoming radiation, but also it can produce masking of the detector by allowing 
exposure of only a certain areas of it. Note that we ignore the effects of the wavefront 
tilt and uneven illumination of different pixels on the detector, and defer the discussion 
of these issues to a subsequent publication. 

Assumptions above result in the following expression for the density of the photon- 
counts registered by the detector: 



where we corrected the total instrumental throughput for the quantum efficiency of 
the detector, J-'(k) = a(k)J-'o{k)- (This form of the notation is quite sufficient for the 
error propagation and sensitivity analyses that we will report elsewhere.) Therefore, 




(24) 
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we derived a model for the interferometric fringe pattern that accounts for a number of 
effects of light propagating through the instrument. In the next Section we will derive 
observational equation to be used for estimating the fringe parameters if interest. 

3. Parameterization of a Polychromatic Fringe Pattern 

As we see from the previous discussion, description of the interferometric pattern in 
the polychromatic case that involves a finite bandwidth of radiation - is a technically 
complicated task. Thus, the observational conditions in the case of polychromatic light 
are significantly altered compare to the simplicity of the monochromatic process. In 
general, all the quantities involved are complicated functions of the wavelength. In the 
previous Sections, we choose a way to describe this process is to collect contributions 
of all infinitesimal constituents of polychromatic light at different wavelengths within 
the bandwidth of the incoming electromagnetic radiationP 

The total number of photo-electron counts, N, registered by a CCD detector per 
wavenumber and per unit time, may be given by the following expression: 

dN(k,t) = F(k)l (k)(l + V(k) sin [0(Jfe) + kxty^dkdt, (25) 

where F{k) is a dimensionless factor representing the total instrumental through- 
put J^^ 1 To(k), V(k) and <p(k) are the spectral density of photon flux, visibility and 
phase of the incoming light; x(t) is modulated internal delay. We are using a nomen- 
clature where a wavenumber k relates to the wavelength as follows k = We also 
accounted for the nominal ~ phase shift due to the SIM beam splitter, which produces 
a sine fringe rather than a cosine one. 
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Note that the total instrumental throughput depends on a number of other fac- 
tors, some of these are the collective area of the detector, quantum efficiency of CCD, 
and overall spectral response of the instrument (as discussed in Section HE} . Our goal 
here is to derive observational equation that may be used to estimate the apparent 
fringe phase and visibility. To estimate the true source visibility and phase one would 
have to perform a set of additional calibration and estimation procedures that will 
be addressed elsewhere. 

A. Integration Over the Spectral Bandwidth 

In this Section we will perform integrations of Eq. (|23|) over wavenumber space and 
time, that are necessary to derive analytical model. This model will be used further 
for the purposes of the fringe parameters estimation. 

Let us first perform integration over the SIM wavenumber bandwidth k G [A^ M , kg m ], 
where fc^ IM = 450 nm is the beginning of the SIM bandwidth, and k$ m = 950 nm is 
the end of this bandwidth, thus k G [450,950] nm. A formal integration of Eq. ([25|) 
over dk leads to the following result 

N(t)Ak sm = / N(k,t)dk, Ak sm = k£ m -k sm . (26) 

^ ^SIM 

In the case of channeled (or dispersed) spectrum output, the integration of this equa- 
tion over the range of wavenumbers is straightforward. For this purpose, we designate 
index, £, to denote a particular spectral channel. Suppose that there exists a total of 
L spectral channels, thus £ G [1, L\. 

Our definition for the spectral channel i implies the width of the channel Ake = 
h\ — kj and existence of a "central" wavenumber ke within this channel. We also, 
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assume continuous spectrum within the bandwidth, so that there is no gaps exist in 
the interval k G [ks m , k^ m ]. A consequence of this is the equality kj +l = k^, which 
leads to the following discrete representation of the bandwidth AA; S im = J2i=i Afcf. 
These assumptions allow us to rewrite the right-hand side of Eq. (j2T?j) as follows: 

/ N(k,t)dk = J2 N(k,t)dk = J2 N t (t)Ak t , (27) 
J k sm i=\ e=i 

where Ne(t) is the instantaneous number of photons within a particular spectral 
channel and has the following form 



N e (t) = f " F{k)l {k)(l + V{k) sin [<P(k) + kx(t)])dk. (28) 

This equation is our first important result. It will help to focus our attention from the 
discussion of coherent processes within the whole wide bandwidth, onto addressing 
this processes on a smaller scale - within a particular narrow spectral channel, t. 

B. Definitions for the Fringe Parameters 

At this point, it is convenient to introduce a set of useful notations. First of all, we 
define the average total intensity of incoming electromagnetic radiation, I Q i, within 
the £-th spectral channel as 



= TT / F(k)X (k)dk. (29) 
Afc* Jk~ 

It is natural to introduce normalized intensity of light Xqi within the £-th channel: 

Mk) = ^ k)Uk) with -L f"%(k)dk = i. (so) 

These new notations allow to present Eq. (|2*K|) as given below 

N e (t) = I °e( 1 + ^f k l %i{k)V{k) sin [<j>{k) + kx(t)] dk^j . (31) 
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In the next Section we will define the fringe visibility, phase and mean wavenumber. 
C. Fringe Visibility, Mean Wavenumber and Phase 

To further simplify the obtained equation, we will introduce functional form the fringe 
visibility, the phase and the wavenumber notations. Thus, the fringe visibility, V^, 
within the £-th channel is given as 



V M = aTT / lot(k)V(k)dk. (32) 
Afc/ J k - 

Similarly to Eq. fl30[) we denote normalized visibility in the channel as 

y %t{k)V{k) _ Hk)Uk)V(k) 1 [% w . 1 

Voe(k) = = , — / V w {k)dk = 1. (33) 

These definitions help us to re-write equation (j3*Tj) in the following compact form 

N t (t) = loe (i + Vo^J^ Mk) sin [<f>(k) + kx(t)] dk) . (34) 

To define mean wavenumber, ke, and mean phase, <f>£, for the £-th spectral channel 
we will use the following expression: 

k e = -L [ ke i oe (k)kdk = -J— [ ke F(k)l (k)kdk, (35) 
&k e J k - loe^ki J k - 

There are two ways to define the phase within the channel. Thus, it is tempting to 

define the mean phase as 

fo = Tu I e ioe(k)<P(k)dk = —j— \ " T{k)Z^(k)^(k)dk. (36) 
Ak< J k - ToAh J k - 

This definition is acceptable for narrow spectral channel, however for a wide channel 

one needs a more convenient form, namely 

4>{ke), which is <f)(k e ) ^ <p e , (37) 
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and is simply the phase value at the central wavenumber. In our further analysis we 
will be using this latter definition. (The relationships between the two definitions for 
the phase Eqs. (jHHjl and (pT7jl will be addressed in Appendix IXj) . 

The three introduced quantities (i.e. visibility, mean wavenumber k# and phase 
at the mean wavenumber (j>{ke)) allow to proceed with integration of Eq. (JHH)- 

D. Complex Fringe Envelope Function 

Definitions introduced in the previous Section allow us to separate functions ki and 
4>(ke) from the functions with direct dependency on the wavenumber k. As a result, 
Eq. (jUIJ), for the total photon count, may be presented as below 



N e (t) = l Al + 



1 



k , 



+ % sin [<j)(k e ) + k e x(t)} — j V oe (k) cos [{k - k e )x{t) + <f>(k) - </>(k t )]dk + 



1 



+ V oe cos [<t>(k e ) + hx(t)] — V 0£ (k) sin [(k - k e )x(t) + <f>(k) - 4>(k e )] dkj . (38) 

To further simplify the analysis, it is convenient to introduce the complex fringe enve- 
lope function, W^[Afc^, 4>{kz), x(t)] , which is given as 

Wl[Akt,(f>(ki),x(t)] = ^£> ,(A;)e J (( fc -^W + ^-^))d^ (39) 



h 



F{k)l Q {k)V{k) e j ((k-h)x(t)+m-m)) dk _ ( 40 ) 



As a complex function, We (please refer to discussion of the fringe envelope function 
given in Appendix [Bj) niay be equivalently presented by its real, Re{lT^}, and imaginary, 
lm{H^}, components: 

We = Re{W t } + j \m{W e } (41) 
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with 



1 



Re{W e } = — ^ Voe(k) cos [(k - k e )x(t) + <f>(k) - <f>(k e )]dk, 



1 



N e {t) = l oi (l + V m sin [<f>(k e ) + A*a?(t)] Re{w< [*(«)] } 

V oe co S [(ftkt) + k e x(t)] \ m {w t [x(t)]}\ (43) 



lm i^} = XTT / Voi(k) sm[(k-k e )x(t) + <f>(k)-<j>(k e )]dk. (42) 

This definition of complex envelope function given by Ea.(|39 )) -(|42 )) allows us to present 
expression (|38|) in a simpler form: 

^sin [<f>{ki) + k e x(t)] Re{w t [x(t)] } + 

The complex fringe envelope function, W^JWi)], as any complex function, may also be 
represented by its amplitude and its phase, namely: 

W t [Ak e ,<f>(k t ),x(t)] = £ e [Ak e ,^k e ),x(t)]e^ Ak ^ k ^l (44) 

where Eg and 0,£ are the amplitude and phase correspondingly. For the complex envelope 
function Eq. (|41|) these two are given as follows: 

S t (t) = ^Re 2 {W £ } + \ m 2 {W £ }, Qe(t) =ArcTan{ ^||^| }. (45) 

Finally, we re-write Eq. (|43[) in the following general form: 

N t (t) = l oe (l + V oe £ e (t) sin [<f>(k e ) + k t x(t) + n e (t)] ^ . (46) 

Note that the apparent visibility of the fringe now is the product of the true averaged 
visibility and the modulus of the Fourier transform of the filter function, evaluated at the 
current delay or 

f x = ViuWi = V i E e e 1 . (47) 
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It is known that the transfer function Wt describes the coherence envelope.^ If Voe(k) ~ 
J-(k)Io(k)V(k) is symmetric, then Wg is real valued, = 0, and only at zero delayf^where 
the envelope is at peak, is the true visibility observed. 

E. Temporal Integration 

The last integration to be performed in Eq. (j25[) (or equivalently Eq. (|4(i|) ). is the integration 
over time. The optical pathlength difference may be modulated either as a set of discrete 
values corresponding to a number of steps in the OPD space (stepping OPD modulation) 
or by ramping OPD over the range of values. The total integration time, At, is the sum 
of durations of eight temporal bins. (While our result is applicable for arbitrary number of 
temporal bins, the SIM design will utilize 8 temporal bins): 

Af=8 

At = t + - 1~ = ^2 An, with An = tf-tr. (48) 
i=i 

Direct integration of Eq. (|43|) leads to expression for the total number of photons collected 
at each stroke of OPD modulation: 

t + N=8 r t+ N=8 

N e At = / N e (t)dt = / l N e (t)dt= Y.NaAn, (49) 

i=l A i=i 

where N^Ati is the total number of photons collected in a particular i-th temporal bin and 
for the £-th spectral channel. Substituting Ng(t) from Eq. (|46jl directly into Eq. (|49|). one 
obtains following expression for 



1 Hi 

Na = ^ ^ N t (t)dt = 

+ V oe sm(f>(ki)£i(t) cos [k e x(t) + £l e (t)] 
+ V i cos 4>{ki) £e(t) sin [k e x(t) + n e (t)] ) dt. (50) 



i r* 



At, Jt 
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To complete this integration, we assume that quantities X e, Voi, (f>(ke), and ki do not 
change with time during the photon-counting intervals. The only quantity that is 
explicitly varies with time - is the optical pathlength difference x(t). 

What makes the polychromatic case much more difficult to study is that the 
different constituents of the finite bandwidth are not responding to the phase modu- 
lation in the same way due to a short coherence length, as shown in Figures HI and [7| 
In particular, Figure |H1 shows a typical interference pattern of the polychromatic light 
that was composed as a superposition of a number of monochromatic constituents. 
Also, in the upper plot of Figure [7| we show behavior of three monochromatic sinu- 
soidal fringes plotted as functions of OPD. The OPD is modulated in equal steps over 
a wavelength as shown in the upper plot of Figure EJ Note that largest wavenumber 
(plotted by thin dashed line) produces the faster changing fringe amplitude. Moreover, 
the polychromatic fringe composed from these three waves does not exactly repeats 
the behavior of either of its constituents. As a result, one would need to minimize the 
bandwidth Ak in order to be able to describe the polychromatic phenomena. 

The integration over time may be performed in a general form and corresponding 
expression for the photon count, N a , is given as follows: 

N u = 2o/l +V oe sm<f>(k e )Re{V ei } + V oe c^s<f>(k e )\m{Pe i } > j, (51) 

with quantities Re{P&} and Im!^} given as below: 

i r*t 

Re{Vu} = — I £ t [Ak e , </>(k e ), x(t)] cos [k e x(t) + n t [Ak t ,(/>{kt),x(t)]]dt, 

i r*t 

lm {A} = 7— I £i[Ak e} <j)(k e ) } x(t)]sm[k e x(t) + tt e [Ake,<j)(k e ),x(t)}]dt. (52) 
At, J t - 
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Optical pathlength difference, x 



Optical pathlength difference, x 



Fig. 6. Shown on the left are three monochromatic waves within the narrow 
finite bandwidth. Right plot shows interference of the polychromatic light that 
was composed as a superposition of these three waves. Note the envelope cor- 
rection to the fringe visibility as a function of the OPD. 

For convenience of further analysis, we combined these two real-valued matrices into 
one complex matrix Vu 



Re{V £l } +j \m{V a }. 



(53) 



Furthermore, with definitions for ReiP^j} and Imj^i}, Eq. (J52)) . this complex matrix 



may be presented as given below 



1 




which, with the help of Eq. (|44p. is conveniently transforms as follows 




(55) 
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a 




o 



o 



Optical pathlength difference, [bins] 




Optical pathlength difference, [bins] 



Fig. 7. Upper plot shows typical behavior of three monochromatic sinusoidal 
fringes modulated in 8 equal steps over a wavelength (as shown in upper plot 
of Figure EJ) . Parameters for the waves are given as <fi = kiX, 0o = 0, with 
k\ = 0.75&0 (thick dashed line), k 2 = 0.8/co (thin dashed line), k 3 = 0.85k 
(thin solid line), thus the width of the spectral channel is Ak = 0.1k . Lower 
plot shows polychromatic fringe composed from these three waves. 
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where the complex envelope function We given by Eqs. (|35|) - (|30|l . 

Eq. (}55|) may further be transformed to establish its true dependency of the 
integrand on time and wavenumber. To do this, we substitute expression for complex 
envelope function, Eq. (|3Uj) . directly into Eq. and obtain the following form for 
matrix Vm with explicit dependency of the integrand: 



1 



F{k)l {k)V{k) e j (**(*)+*(*)-*(*<)) 



dtdk. (56) 



y Nin j 



1; ImjPfl}; Re{Va] 



J 



loeVoiAktAri 

Finally, we have defined everything that is needed to study Eq. (J5T|) . for the 
polychromatic fringe, which may equivalently be presented in a matrix form as below 

<f>(k e ) I ' ( 57 ) 
1; \m{Vw}] Re{V iN } J y IqiVqi sin <j>{ki) 
where the complex matrix Va is given by Eq. (J56|) . 

The obtained result Eq. ()5ip. (|5bp (or, equivalently Eq. (|57)l) constitutes the gen- 
eral form of expression for the polychromatic fringe. We will use this result to finalize 
the development of the general from of the observational model for polychromatic 
case with arbitrary phase modulation. 

Ideally, one would need to determine not only three quantities X ^Vo^ cos 4>(h), 
JoiVoi sm(ki) , and T oi , but the full functional dependency of the original quantities. 
However, the finite width of the observational band-width Afc^ complicates the es- 
timation process by bringing the non-linearity in the observational equation via the 
envelop function W. Note that if one neglects the size of the bandwidth Ake with 
respect to the mean wavenumber kg or Ake/ke — > (the envelop function becomes 
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unity W — > 1), one recovers the full simplicity of the monochromatic case represented 
by Eq. ©P 

4. General Solution for Polychromatic Phasors With Noisy Data 

Currently in use, there are two fringe estimators, one for visibility (the unbiased 
estimator is V 2 ), and one for the phase (the unbiased estimator is the complex phasor). 
The V 2 estimator is already worked out in much detail (i.e. Refs. l4J10llT3j) if the 
complex phasor estimator is completed. So the development of the complex phasor 
was the main purpose for the presented work. As it is known, the complex fringe 
visibility can be represented by a phasor; if the fringe is stable, we can add the phasors 
vectorially over multiple samples. This co-adding can provide an improved signal-to- 
noise ratio. To co-add the fringe phasors requires a phase reference, for instance the 
white light phase !^ 

In this Section we will develop optimally-weighted solution that accounts for a 
number of noise sources and will be applicable for a general case of delay modulation. 

A. Fringe Equation for Noisy Data 

For the purposes of clarity we will omit spectral index i. All the obtained results are 
valid for any channel and thus could be easily reconstructed, if needed. 

In the case of noisy data, observations of photon-counts Ni are actually done 
with errors and, in reality, we observe Ni — Ni + e*, where Ni is the value of photon 
counts at the i-th temporal bin in the absence of noise and is a random variable, 
representing the noise contribution to the measurements. We assume that £j are ran- 
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dom variables that are primarily due to gaussian statistics. (This approach may be 
extended to incorporate other sources of noised The corresponding results will be 
reported elsewhere.) that are distributed around zero and following relations are valid 



Ni = Ni + e< 



EUi) = 0, 



E{4) = °l 



(5? 



In the general case one must account not only for the Gaussian statistics of read-out 
process, but also for the Poisson statistic that governs photo-emission (and photo- 
counting) process. Thus, a correct approach would be to assume that is a sum of two 
terms = /i ef+ (1 — fi) ef, where ef is Gaussian and ef is Poissonian variables and \i 
is a number between and 1. Expected complication arises from the fact that standard 
deviation computed for the photon-counting Poissonian bias is actually proportional 
to the signal (erf) 2 oc Zj (see discussion in Ref. ITUJ) . This issue is out of scope of the 
present paper and we will address this issue at a later time. 

We also assume that N{ are independent, therefore, we may form a diagonal 
covariance matrix for the quantities Ni (or equivalently for e^) with dispersions of on 
the diagonal: 
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(59) 
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where G y is the matrix of weights. Therefore, in the case when noise is present in the 
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data, equation (joTj) has following matrix form as below 



Ni + et = (l; lm{P«}; Re{f> H } 



1 Z ^ 



T oi Voe cos <j>(k e ) 
XoiVoe sin 4>(k e ) 



(60) 



or, equivalently, 

Nt + ei = (61) 

with indexes z and a running as t G {l,...,iV} and a G {1,2,3}. Vector X a = 
^X £; X ^Vo£ cos0(/^); ToeVoe sin <K^) j is the vector to be determined and matrix 
A T = Ai a = (l; ImjP^}; Re{P&}) is the 3xiV rotational matrix in the phase space. 
A maximum likelihood solution to the system of equations ()61|) may be given by the 
following system of equations 



N 



X a = J2 A o a Ni, where = (A T G y A)~ 1 A T G y , (62) 

i 

with Al being an optimally-weighted pseudo-inverse matrix. Note that by choosing 
different gain matrbP instead of optimally weighted least-squared matrix Eq. ([59)1 . 
one may obtain solution with different, specifically designed properties. Nevertheless, 
our solution has enough embedded generality as it allows for arbitrary properties of 
noise contribution, which will be further explored below. 

B. Optimally- Weighted Pseudo-Inverse Matrix 

In this Section we will find solution for the pseudo-inverse matrix Aj, that was intro- 
duced by Eq. (j62|) . To construct this matrix we will use the weights matrix G y given 
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by Eq.()59p and the rotation matrix A given by Eq.(J33J) as: 



A, = (l; ImjP,}; RejP,}) = (l; s f ; q), 



(63) 



where we denoted Si = lm{Pi},Cj = Re{"Pj}. 

Let us construct matrix (A T G y A) first. Calculation of (A T G y A) is straightfor- 
ward even for the most general case of arbitrary number of temporal bins (N > 3) and 
with arbitrary integration intervals (Ari 7^ Atj for i ^ j). Thus, after some algebra 
we find the following structure: 
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(64) 



i 1 i 1 i 1 

By inverting the obtained result one constructs the covariance matrix A of the fol- 



/ 



lowing structure: 
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where determinant of the matrix (A G y A), A = det ||(A G y A)||, is given as 



N 



A, 



-Y 

ijk 



y S{Cj SjCj,j 



(65) 



with a triple summation for all the indexes denoting the temporal bins and running 
from 1 to N, namely V k} G [1, ...,iV]. 
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These intermediate results allow us to write the solution for the (iVx3) optimally- 
weighted pseudo-inverse matrix Aj, = A\Z in the following compact form: 



At = (A T G y A)- 1 A T G y 



(66) 



where coefficients A%, B%, C£ and V° depend on duration of each temporal bin, mean 
wavenumber and variances for the data taken in each bin, and are given by 

N 1 

A k ^ 2 2 2 {j^i^i ^j^i) SjCi) ~\~ (SjC k -SfcCj) -|- (s^Cj S{C k 

ij a i a j a k 

N 1 

®k = 2 2 2 ( C » _ C i) [( S i C i ~~ S i C *) + ( S 3 C k ~ S k C j) + ( S k C i ~ S i C k)] , 

ij i J ft 

N 1 

C k = ~ Yl ^2T2T2 ( S i ~ S i) l( S i C 3 ~ S 3^) + ( S 3 C k - S kCj) + (s k 



N 1 N 1 2 

X>° = ^ ^ = 77 ^ ? ? s [( s i c j - s j c i) + ( s j c k ~ S k Cj) + (s k Ci - SiC k )] . (67) 
k ijk * 3 k 

Definitions for the quantities Sj and Cj, Si = Imj'Pj}, Ci = Re {'Pi}, allow to present 
expressions (16 7|) in a more convenient form. First, remember that complex matrix, Pi, 
as any complex function, may be represented by its amplitude and its phase, namely 



Pi = Re{Vi}+j \m{Vi}= Pie?** 



(68) 



where Pi and 7Tj are the amplitude and the phase of this complex matrix correspond- 



ingly and are given as follows: 



Pl = ^Re^J + lm 2 !^}, vr, = ArcTan{^rj> 



(69) 
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with complex matrix Vu is given by Eq. (|55j). These quantities allow presentation of 
Eqs. (167)1 in the following form: 



At = E 



at . r , 

PiPj Sin[7Tj - 7TjJ 



x sinfvri - 7Tj] + sin [7^ - n k ] + p k p t sin[7r fe - 7Tj]J . 

K = E 



AT / \ 
(Pi COS 7Tj — COS 7Tj 1 



x fp^- sin[7Ti - 7Tj] + pjp fc sin^ - vr fc ] + p k p { sin[7r fc - 7T;] Y 

N 1 ■ \ 
[Pi Sin 7Tj — Sin 7Tj J 



/ / % j k 



X 



A' 



x (piPj sin [7^ - 7Tj] + Pjp k sin [7^ - vr fc ] + p fe p; sin[7r fc - 7Tj] J , 

1^1 2 
= 2 2 2 (PiPj _ ^i] + PjPk sm[7Tj - ix k ] + p k pi sin[7r fc - 7r»] ) . (70) 

At this point we have all the expressions necessary to present the optimally- 
weighted solution for the polychromatic phasors. 

C. Photon Noise- Optimized Solution for Polychromatic Phasors 

An optimally-weighted solution for the quantities X a may be obtained directly now 
from Eq. (|b2|) with the help of expressions (|66j) - (|67j) in the following compact form: 



1 N 



N 



k 
N 



l« V °sm4>« = ±.J2N k Cl (71) 
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with coefficients of B%, C£ and V are given by Eqs.()67p and (|70j) . 

The obtained solution for the polychromatic visibility phasors given by Eq. (|7T|) is 
given in the form of a linear combination of weighted photon counts recorded during a 
particular integration period. This form turned out to be very helpful when analyzing 
contributions of CCD pixels that are systematically biased. The obtained result may 
be used to de-weight 'bad' pixels (in a statistical sense) and, thus, to reduce the 
problem of biases while estimating fringe parameters. 

This form allows to express an optimally-weighted solution for visibility, phase 
and the constant intensity terms in a familiar compact form: 



N „ N _ 9 

{12 Nk&l) + (J2 N k CQ 




k k 



The form of the obtained solution is simple to understand and it is straightforward 
to implement in the software codes. All the information necessary to calculate the 
3iV coefficients of A%, B%, C£ and V is presumed to be known before the experiment. 
Thus, for the case when N = 8 one would have to calculate only 24 numbers from Eq. 
(JHTJ) . These numbers correspond to 8 numbers of A%, 8 numbers of £>£ and 8 numbers 
of C£. The experimental data is used as input to Eqs. (fTTj) (or directly Eqs. (JZ2J) to 
produce the best estimates of the actual values for visibility, phase and the constant 
intensity term. This approach is currently being utilized and corresponding results 
will be reported elsewhere. 
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5. Filtered Polychromatic Light: Spectral Channels with Narrow Bands 

It is thought now that SIM will be operating at 80 spectral channels for the science 
interferometer and at 4 channels for both guide interferometers. The data will be read 
at a millisecond time rate. SIM will be dispersing light just before the interfering 
radiation reaches the detector. To take this fact into account, we define a filter that 
allows to limit the total bandwidth of the incoming radiation. We will formally denote 
a filter with such a properties as follows: 

Ti{k) = Ti{k - he; Ak e ). (73) 

We assumed that this filter operates within the £-th spectral channel by allowing to 
pass through only such a radiation that is composed from the frequencies correspond- 
ing to the mean (or central) wavenumber k e within the bandwidth of Ak e . Thus, filter 
J-'e(k) enables the instrument to "see" light only in a certain interval Ake around the 
mean number ki. 

In this Section we will develop a model that employes such an approach for the 
case of a rectangular bandpass filter. Due to its analytical simplicity, the rectangular 
bandpass filter is the most known construction in the Fourier optics. This analysis 
will allow us to establish correspondence with the previously obtained results both 
for monochromatic and polychromatic light. 

A. Rectangular Bandpass Filter 

To take advantage of the results derived in the previous section, we must first decide on 
the properties of the bandpass filter. This decision in return will affect the properties 
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of the envelop function. Below we shall develop a model for a special case of the 
bandpass filter - a rectangular bandpass filter denoted here as Tt, which is done 
analytically in the following form 

Hk) = $>(*), where F t {k) =\ = ^ * G [ *'~ Kl (74) 

[ o, kt[kj,kt\- 

We can also assume that the width of a spectral channel is small, so that both intensity 
of incoming radiation, X (fc), and apparent visibility, V(k), do not change within 
the spectral channel (in particular, this leads to Vae(k) = 1 in Eqs. and 
Therefore, the following conditions are satisfied with a particular spectral channel, £: 

Xq(/c) = const, = const, JF ^ = const, (75) 



<j>{k) - <f>{k t ) = dot(k-k e )+0(^), (76) 



where doi = is the delay within the £-th channel. 

One may perform integration of the fringe envelope function which is 

given by Eq. (or use equation (|112|) for the unperturbed envelope function and 
then apply iterative procedure outlined in Appendix[Bl). To the second order in phase 
variation (i.e. C(|^|)), the resulted envelope function has following properties: 

W e [Ak e ,<p(h),x t ] = - 1 /*V(fc)2b(fc)V(A ; )^&*-* < ^ + *«-^)dfc = 



ZofV^Afc^ j k - 



Afe> 7fc- dkj ■ 



Let us introduce a convenient variable, k = k — kg, and remember that Afc^ = 
fe/ — kj and = |(A;/ + fej") are the width of the spectral channel and the mean 
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wavenumber. This allows us to integrate expression (|77jl over the wavenumber space 

w t[ A k „ M - w F^^- ^VT - (78) 

AkiJ_i Ake ^Ak e (doe + x(t)) 

We can now present Eq. (J55j) for matrix Vu as follows: 

1 m sin[±AA>(GL> + x(t))} 

At this moment, we introduce another convenient variable, r = t — ij. Analogously, 
At; = t+ — £7" and ^ = + tj 7 ) are the duration of the temporal integration within 
the i-th bin and the mean time for this bin correspondingly. This result is used to 
transform Eq. (|79|l as below: 

Pa = e jk ^ SV a , (80) 
with matrix coefficient SVu given by 

AtJ_1 At( \t\k,(d <X + x(t i + T)) ^ 

The obtained result explicitly depends on the functional form of the OPD mod- 
ulation, x(t). To integrate this equation one first needs to make certain assumptions 
on the temporal behavior of x(t), which will be done in the following Sections. At 
this moment we present Eq. (|51J1 in the following form: 



N u = Zoi\l + Votwn [<f>(ke) + k t x(U)] ^e{SV B } + 

+ V oe cos [<f>(kt) + hx(U)] \m{SVu}j , (82) 

where the complex matrix SVu is given by Eq. (|81)). The importance of separat- 
ing terms with SVu is i n the fact that one can establish clear correspondence with 
monochromatic light, for which SVu — I Hi the identity matrix. 
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In the next two subsections we will study two different special cases of OPD mod- 
ulation, namely the stepping and ramping modulations of the optical path difference. 

B. Stepping Phase Modulation 

The stepping phase modulation realized when the pathlength difference is changes 
as a set of discrete values corresponding to a number of steps in the OPD space, 
mathematically this process represented as follows: 

x(t) = x(tj), where x(U) — < " I ■ ' ! (83) 

| o, t t[tT,t+}, 

with ti = \{tf + tj 7 ). This procedure defines the temporal bins that will be used to 
modulate the interferometric pattern. 

Conditions (}83|) allow for a significant simplification of temporal integration in 
Eq. (fRTjl . It simply is leading to a substitution x(t) — > Xi in Eq. (JEJ), and matrix Vu 
takes the following form 

V fA = e ]kex{u) — ^ — r^, (84) 

±Akt(dot + Xi) 



As a result, to the second order in the phase variation (i.e. <p{k() m <pi + 0(-^ )), 



observational equation (J82|) is taking form as below: 



N u = 2 oe ^1 + Vot sinc[~Afo (x; + doe)] sin [</>(h] + hxi] ^ . 



55) 



The obtained result Eq. ()85|) clearly depends on the particular form of the envelope 
function. As such, it has most of the parameters that are necessary for the phase 
estimation purposes in the case of wide bandwidth. 
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For the most practical cases the value of the sine function will be close to sine ~ 1. 
Indeed, let us analyze the argument of this function, ~Ake(xi + d e). Thus, one might 
expect that within the spectral channel the phase will stay constant, hence doe = 
J^j ~ 0. Furthermore, for the estimation purposes let us assume that all the step- 
sizes %i are essentially Xi — i^, where N is the total number of temporal integration 
bins, i is the number of a particular temporal bin, i e [1, N], and Ao is the reference 
wavelength chosen for the OPD modulation (or Ao = f^-, where ko is the wavenumber 
corresponding to chosen modulation wavelength). Also remember that the width of a 
spectral channel is related to the total SIM bandwidth as Ake = A ^ IM , where Ak$\M 
is the total SIM bandwidth and L is the total number of spectral channels used for 
the white light fringe detection. Therefore, one has 

1*1/ 1 a Ak$\M 

-Ak t {Xi + d oe ) rs -Ak £ Xi = -7-77—; • (86) 

Z Z Lily Kq 

Assuming Ag| M = 450 nm and Ag !M = 900 nm, and Ao = 900 nm, thus yielding 
A ^ IM = 1. The maximal value for the expression (|86|) is realized when % = N, thus 

Tri AA; SIM ^ 



LN k ~ L 

Currently, there are different numbers of spectral channels used to process data from 
our testbeds. This number may be as large as L = 80 and as small as L = 4. Of coarse, 
when L = 80, the ratio n/L becomes n/L = 0.03927 and, thus, s\nc[^AkeXi\\L=80 = 
0.99974, and similarly for L = 4 the sine function becomes sinc[|AA;£Xj]|L = 4 = 0.90032. 
We will address the issue of phase estimation sensitivity to the width of a spectral 
channel Ake at a later time. 
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This observation allows us to present the sine function as a series with respect to 
the small parameter Ak(,z (with z being defined as z = x-i + doe) as 

A/c £ Z ^ (-1)™ rAkiZnZn 



n=0 



(2n 

3 , L 2 J 5! L 2 J l 7! L 2 J J 



one can present the expression Eq. (|53jl in the following form: 



1 + 



+ Vb, (1 ^4 ^ + £V 192Q '-) sin + keXi] j . (89) 

The obtained expression models the expected number of photons detected at the 
CCD for the rectangular bandpass filter and stepping phase modulation. It extends 
the results obtained for the monochromatic case on the finite size spectral bandwidth. 
This fact is indicated by the explicit dependency of the obtained result on the width 
of a spectral channel Ak£. (For the most of the interesting practical applications, the 
size of the delay within a particular spectral channel is very small doe = §^ ~ 0, 
which further simplifies Eq. fl89j) ). 

C. Ramping Phase Modulation 

In this Section we will discuss another type of phase modulation - the case when 
the phase is linearly changes with time. This modulation utilizes the phase ramping 
technique. (For more details, see RefsEE) To develop analytical solution we will be 
using the system equations developed above, specifically Eqs. (JoTj) and (|55|). 

The optical path difference for the case of ramping phase modulation is modeled 
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as a continuous function of time as follows: 



x(t) = x + v t, 



(90) 



where xq is the initial OPD value and v is the instantaneous velocity of OPD mod- 
ulation. Remembering the definition for r as r = t — ti, and Afj = tf — t~[ and 
ti = \{tf + tr), Eq. (jHU) takes the form: 



(91) 



with coefficient SVa given by 



1 r + ^ An •, sin[|AWr)] 



(92) 



and z(t) = + x(tj) +vt. This allows us to present Eq. (|82j) in the following form: 



N u = Xoi[l + V oe sin [</>(k t ) + k e x(U)] Re{5T ei } + 
+ V^cos [(j>(kt) + k e x(ti)] \m{5Vei} ), 



(93) 

where the complex matrix of additional rotation in the phase space, 5V&, is given by 
Eq. (J81|) . Equation (J93j) may equivalently be presented in a matrix form as below 



N, 



1; sin k£x(ti); cos k(x(ti) J x 
/ 



x 



. (94) 



1 

Re{5Pn} -\m{5P u } 

lm{5P«}; Re{5V ei } J \ J, oe v tsm(p{Ke) J 
The result of integration of Eq.([92jl may not be presented in a compact analytical 
form. It rather could be expressed in the form of two functions defined as Sinlntegral 



loiV oe cos </>(kt 
loiVoe sm<j)(k e 



\ 
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and Coslntegral. To simplify the analysis, the sine function may be given in the form 



of power series expansion with respect to the small parameter Aktz(r) as given by 
Eq. (|HHJ) - This expansion allows us to present Eqs. (f9*2"|) in the following form: 



SVm 



1 f^ An ., sin[iAfc £ z(r)l 



0(^)1 dr, (95) 



21 1920 ' v 7!2 6 

where z(r) = doe + xfy) + vr = Zi + vr with Zi = d e + x{ti). This equation, (|9*5jl 
was integrated to obtain the following result for 8Vu: 



- sin[^Ar t ] - Akj ~ Ak\ Ak* 

1 Hi — ~ 1 r A.U 77779 + °£i '-, , a + 



\kov At,- 



24ft 



1920 kt 



•7! 2 6 A^ 



(96) 



where complex coefficients Au is given as follows: 



A. 



l + jhzi)) 2 -(~k £ v An) 2 



- 2 (l - jkgz^j cos[iA; £ u At;] 



sin[|A;^ v At^ 
\k t v At { 



(97) 



and £>£j is computed as follows: 



B 



n 



1+ (l-jkezi) 2 -(-kevAn) 2 



+ A{2 - jk iZi ) 2 + 



1 v „\ \ sm\\kev At, 

'2 / / ^kivAr, 



4^5 + (1 - jfe)((l - jfc^i) 2 - (i^f Ar,) 2 ) jcos[i 



ft f Ar,- 



(98) 



The obtained expressions may be used to simplify the results of temporal integration 
Eq. (|9*3|). As a result, the coefficients Re{5P#} and I m { (5^^ } in the fringe equation 
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Eq. ()93|) may be written in the following form: 

Re{5P u } = if t J 1 + (2 - $z? - (-k t v At, 



24fc 



2 C o4M>Ar,]|§ + (^k (99) 



lm{ ^,J =%(-^ + 4.A^ +0 (i), (100) 

with Zi = doe + x(ti) = <i £ + x + f £«• 

The obtained expression models the photon flux detected at the CCD for the case 
of rectangular bandpass filter and ramping phase modulation. It extends the results 
obtained for the monochromatic case on the finite size of spectral bandwidth. This 
fact is indicated by the explicit dependency of the obtained result on the width of a 
spectral channel Akg. 

6. Discussion and Future Plans 

The main objective of this paper has been to introduce the reader to the concepts and 
the instrumental logic of the SIM astrometric observations, especially as they relate to 
estimation of the white light fringe parameters. The set of formulae described herein 
will serve as the kernel for the future mission analysis and simulations. We have 
also developed a set of expressions that may be used for fringe visibility and phase 
extraction for both SIM science and guide interferometers. The obtained expressions 
depend on the effective operational wavelength of OPD modulation, the width of a 
particular spectral channel Ahe with the mean wavenumber ki and corresponding 
wavelength A^. Our model accounts for a number of instrumental and physical effects 
and is able to compensate for a number of operational regimes. 
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The logic of our method is straightforward: one first assumes the desirable prop- 
erties of the bandpass filter, then finds the corresponding envelope function, and then 
applies the obtained expressions (which are valid for a generic case). The obtained 
solutions for the envelope function W and, most specifically, 5V& may be directly sub- 
stituted either in the expression for the complex visibility phasors Eqs. ()7()jl - ()71j) , or 
into equations for the visibility, amplitude and phase of the fringe, given by Eqs.(|72|). 
We applied this formalism to the case of a rectangular bandpass (the obtained re- 
sults are given by Eqs. (|118|) - (jl21j0 . While the complex visibility phasors are linear 
with respect to photon counts, the explicit expressions for the fringe parameters are 
non-linear. This fact may be used to design specific properties of unbiased fringe 
estimators for processing the white light data. 

Having developed algorithm to obtain the best estimates for the fringe parame- 
ters, Eqs. (I72|) it is naturally to write down the expression that would allow to estimate 
the group delay for the polychromatic case. It turns out that the following expression 
d = jf^2e{4>i/ke) is sufficient to estimate the group delay for the SIM bandwidth 
to the required accuracy of a few tens of picometersP^ This single channel error ex- 
pression can be used to determine the group delay error when combining multiple 
spectral channels of data via phase delay or group delay methods. These errors will 
vary based on the underlying assumptions, but a reasonable figure of merit to keep 
in mind is that a random 1% error in the wavelength combined with the nominal 10 
nm rms delay requirement leads to approximately 60 pm of delay error when using 
4 spectral channels of data processed with a least squares algorithm for each chan- 
nel. This result scales linearly with the wavelength error and the delay offset. In our 
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further work we will numerically address the problem of unbiased estimators for the 
fringe phase, visibility and group delay. This effort is currently underway. 

Our analysis showJ^' that, while the model of the rectangular bandpass filter 
is working quite well, for the 'real life' one must account for the effect of leakage of 
light. This effect due to the leakage of light onto the studied spectrometer pixel of 
the detector from the adjacent pixels with different wavenumbers. At this moment, 
it seams more appropriate that a combination of the rectangular bandpass filter with 
additional effect of light leakage from the adjacent pixels that must be included into 
the model of a CCD detector. The corresponding analysis, simulation results and 
implications for the instrument design will be reported elsewhere. 

A. Two Definitions for the Fringe Phase 

In this Appendix we will address the issue of the mean phase definition which requires 
some additional work. It is tempting to define the mean phase as 

1 X w (k)<p(k)dk = / ' ' F(k)I (k)(j)(k)dk. (101) 



Akt Jk- loeAke j k - 

However, one needs to relate this expression to the quantity <p(ki), which is the phase 
value at a particular wavenumber. Assuming that phase <f)(k) is a slow varying function 
of k and, as such, it may be expanded in a Taylor series as follows: 

<f>(k) = <j>{h) + ^(k - kt) + \^\ k (k - hf + O(Akf). (102) 
We can now substitute this formula directly in Eq. (jlOlj) . which results in 

4>t = <K^) + ^S t\t [ ke Hk)Uk)(k-k e ) 2 dk + 0(Ak 3 e ). (103) 
2 dk 2 k^oeAkt J k - 
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Or in other words, the phase value 4>(ke) at a particular wavenumber ke is related to 
the mean phase <pi within the spectral channel with width Aki (i.e. Eq. ljlOl)) ) by the 
following expression 

(2) 

where n\ is the dimension-less second-order moment of wavenumber distribution 
within the spectral channel of interest: 

/if = f l Hk)1 Q {k){k - k e ) 2 dk. (105) 

In the general case, when the higher order moments are considered, this expression 
takes the following form: 

p 

<l>(ke) = AA? /.f +0(A^), (106) 



p\ dkP 



(v) 

with moments [if given as follows: 

•k 



/if = 1, 4 1} = 0, $ = -^^1 F(k)l (k)(k-k e ) 2 dk, (107) 

±oe/lK e Jkj 



$ ] = ~ — TTp+T / e Hk)Mk)(k - h) p dk, < |/zf | < 1, Vp. (108) 
l oi Ak p e J k - 

Note that the wavenumbers within a spectral channel may be considered uni- 
formly distributed, thus prompting to use <p(ki) = (f)e + O(Akj). However, the knowl- 
edge of the second moment y\ may be important in combining the fringe solution 
for the whole operational bandwidth. This question will be addressed elsewhere. 

B. Approximation for the Complex Fringe Envelope Function 

Expression for the fringe envelope function Eq. (J39)) . contains terms that are of the 
first and higher orders of phase variation within the £-th spectral channel. We shell 
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separate these terms by expanding phase (j)(k) in the Taylor series around the mean 
wavenumber kg as given by Eq. (jl02j) . This transforms the argument in Eq. (}39|) as: 

[{k - h)x{t) + 4>{k) - <f>{k e )] = (k-k e )[x(t) + d 0e ]+O(Ak 2 ), (109) 



where dn 



dk 



is the group delay within the £-th channel. In the regime of small 



phase variations within the spectral channel Aki 



dk 



= Akid £ <C 1, we can expand 



the exponential argument in the expression Eq. (J3U|) as given below 

expjj [(* - k e ) (x(t) + d M ) + 0(Ak 2 )^ } = 

l+j(k- k e ) dot + 0(Ak 2 )^ expjj (k - fc*)a:(t)}. (110) 

This last expression may be used to re-write the phase-dependent envelope func- 
tion from Eq.(jnHl) as 



Wg [Aki, fa, x(t 



8 1 r +oc 

l + d 0e ^^ + O(Ak 2 )} V m {k)e l{k ~ kt)x{t) dk. (Ill) 



dx(t) 



Defining the unperturbed fringe envelope function (i.e. that is un-affected by the 
phase variations inside the spectral channel) as below 



+oo 



W e [Aki, x{t)] = / Voe(k) e j dk. 



(112) 



We may present expression (jlll|) for envelope function in the following form: 



Wi[Ak t ,(f>t,x(jt)] = W ( [Ak t , x(t)} + d u W' t [Ak t , x{t)] + 0(Akj), (113) 

where superscript ' denotes partial derivative with respect to OPD 5/5x(t). 

At this point we have established the functional dependency of the envelope 
function, but for the immediate purposes we will be using a generic form for this 
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function, Wt [Akt, x(t)] presenting it only by it's amplitude and phase: 

W t [Ak t ,x(t)] = £e[Akt,x(tj\exp{jn e [&kt,x(t)]} =^exp{j^}. (114) 

Similarly to the expression ()113j) we re-write this result in the following form 

Wt[Ak t ,<f>i,x(t)] = ^expjjfi,} + doe (s' e + j^)exp{j'^} + O(Akj). (115) 

At this moment we show the functional form of real and imaginary components of 
the envelope function. Thus, from Eq. ()115j) one immediately has 

Re|w^ [Aki, fa, Xj\ I = £gcosf^ + c£q£ f £^ cos — Efit t sinfM + C(AA;|),(116) 

lm|Wi[AA;^0/,z<] J = Eg sin Sl t + dot\Et sin + Eft' t cos^J + 0{Ak 2 e ). (117) 

The obtained equation exhibits explicit dependence on the phase variation inside the 
spectral channel given by d w = d<p/dk\ k ^. This issue will be addressed elsewhere. 

C. Solution for Rectangular Bandpass and Stepping OPD Modulation 

In consideration of completeness, we present here a general case solution for an 
optimally-weighted visibility phasor for a rectangular bandpass and stepping OPD 
modulation. In the previous Section we obtained this solution in a general case, there- 
fore, the desired solution may be obtained directly with the help of expressions (fTTj) . 
Corresponding optimally-weighted solution may be presented in the form of Eqs. (|71|) 
and (J72j) with coefficients C% and V° depend only on the size of modulation 

steps Xi, mean wavenumber k, width of a spectral channel Ak, variances of the data 
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particular temporal bin. These coefficients are given as follows: 



N 1 



.Akx; 



sine — 



x 



sine sine 



+ sinc[ 



2 

Akx, 



sinc[ — — — ] sm[k(xi — Xj)] x 

J 1 sin[A;(xj — xj)] + 
sine - ;1 | sm[k(xj — x k )] + 



2 J 
.Akxk 



+ sinc[ ^^ ] sinc[ A ^ ] sin[k(x k - Xi)} 



El / Akxi - Akxj - \ 
„ » » sine cos hx* — sine cos kxj x 

crfajai V 2 2 / 



x 



r Akx 



,Akx 



sine 



[— — ] sinc[— — ] sinffcfo - Xj )\ + 



+ smc[ — — — J smc[ — - — J sm[A;(a;j — x k )\ + 
+ sinc[^^] sinc[^p] sin[£;(:r fe - ^)]J , 



„ ^ ^ 1 / r-^X.ktXi-. — y^^}vXj, — \ 

C k = - 2 2 2 (sinc[— — J smkxi - smc[— — ] smkxj) 



x 



Jsincf^p] sinc[^p-] sin[£;(:rj - Xj )} + 
+ smc[ ^ Xj ] sinc[^^] sin[k(xj — x k )] + 
+ sine [ A ^ Xfc ] sine sin[£;(:r fc - Xi )] 



x 



-Akx 



N 



iA i 



3^ o- 2 o- 2 o- 2 

ijk 1 3 k 



.Akxi 



.Akx,-, 



sine 



sine 



r ^^.kx j i 
+ sine sine 



2 J 2 
Akx k 



sin[A;(xj - Xj)] + 



, rAkx k 
+ sine | — - — | sine 



2 

Akxi 



2 J 2 
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sin A; (:c ~ — Xk)\ + 



sin[/c(x fc - x^] 



where, in consideration of brevity, we omitted index I denoting a particular spectral 
channel. 

The obtained result Eqs. ()121j) clearly depends on the sine envelope function and 
thus it has all the information that is necessary for the phase estimation purposes in 
the case of the wide band-pass. Note, that this result assumes that the phase does 
not change inside the spectral channel. Also, this result directly corresponds to the 
result obtained for the monochromatic case. This may be demonstrated by taking the 
limit Ak/k — > 0, which will lead to recovering the familiar form of monochromatic 
fringe with coefficients At, Bt, C\ and V given as follows: 



At 



sm[k(xi — Xj)] + sin[k(xj — Xk)] + sm[k(xk — , (122) 



2 2 2 ( COS ~ cos kXj) X 

sin[A;(a;j — Xj)] + sm[k(xj — Xk)] + sin[fc(x fc — , (123) 



no 



2 2 2 ( S * n ~ s i n kXj) X 

sm[k(xi — Xj)] + sin[k(xj — Xk)] + sxn[k(xk — Xi)] , (124) 



N 



E- 4 



1^1 2 

= 2 2 2 \sm[k(xi - Xj)] + sm[k(xj - x k )] + sin[k(x k - xM .(125) 

^^-r 'vfc qai 

ijk 1 J K 

This form demonstrates that only the terms with i ^ j ^ k ^ i are producing a 
non-zero contributions to the result, while the terms where at least two of the indexes 
are equal (i.e. i = j or % = k or j = k) will vanish from the sum. 
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A. The Optimally Weighted 4-bin (ABCD) Algorithm 

Results obtained in the previous Section are valid for any number of OPD modu- 
lation steps. In order to show their correspondence to well-known formulations we 
will present their particular form for the case with N = 4, equal step sizes and 
wavelength-matched OPD modulation strokes. This will allow us to obtain results 
that are constituting a four-bin algorithm. 

Thus, in the case of 4 equal step sizes all of |A, the OPD is stepping in increments 
of Xi = ji, i G [1, 4] that correspond to the phase changing in steps of <pi = |i By 
calculating the quantities A% B%, C% and V° in Eqs. (|122j) - (jl25|) and substituting them 
into Eq. (|72|) we recover the optimally-weighted form of 4-bin (or ABCD) algorithm: 



Ni [u\ + 2a\ + aj) - N 3 (lc\ + 4 + aj) + (N 2 + N 4 ) L\ - aj) 

loVcosc^o = ^ j± '- 

2 (of + cr| + <T§ + cr| 



NAcjI + 2a\ + of J - N 2 [a\ + of + 2af ) + [N x + N 3 ) (crf - a\ 
Xo l^sin^o = — 



[Nx + N 3 ) (of + al) + (N 2 + N 4 ) (af + a 3 



4 

T = J -^t ^ : r-^- '-■ (126) 

This solution produces following results for the optimally- weighted visibility, V^, and 
the phase, 
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Ni (<4 + 2of + of) - N 3 [2af + of + of) + (JV 2 + JV 4 ) (of - o ; 



+ 



iV 4 (c7? + 2o| + o| 



iV 2 of + of + 2<t| 



{N 1 + N 3 ) (al-al 



(m + N 3 ) [a\ + of) + (JV 2 + JV 4 ) (<rf + a 



b% = ArcTan 



iV 4 (of + 2of + of ) - iV 2 (a? + of + 2of ) + (N t + N 3 ) ( a\ - a\ 



N x ( a\ + 2of + of ) - N 3 (2of + o| + o| ) + (N 2 + N 4 ) (af - of 



(127) 



This is the optimally weighed form of the most well-known fringe estimation 
ABCD algorithm. The absence of extra f phase in the result O is due to the phase- 
stepping approachPED In the case when experimental noise is absent, this result di- 
rectly correspond to the previously obtained solutions for the ABCD algorithm in the 
monochromatic case. 
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